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\ After an elementary derivation of the "time transformation", mapping a count- 

ing process onto a homogeneous Poisson process with rate one, a brief review of 
Ogata's goodness of fit tests is presented and a new test, the "Wiener process test", 
is proposed. This test is based on a straightforward application of Donsker's The- 
orem to the intervals of time transformed counting processes. The finite sample 
properties of the test are studied by Monte Carlo simulations. Performances on 
simulated as well as on real data are presented. It is argued that due to its good 
finite sample properties, the new test is both a simple and a useful complement to 
Ogata's tests. Warnings are moreover given against the use of a single goodness 
of fit test. 
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52 ; 1 Introduction 

Modelling neuronal action potential or "spike" trains as realizations of counting pro- 
cesses / point processes has already a rather long history [18, 19]. When this formal- 
ism is used to analyse actual data two related problems have to be addressed: 
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The conditional intensity of the train must be estimated. 



• Goodness of fit tests have to be developed. 

Attempts at modelling directly the conditional intensity of the train have appeared in 
the mid 80s [10,2,5] and, following a 15 years long "eclipse", the subject has recently 
exhibited a renewed popularity [11,23, 17,22]. The introduction of theoretically mo- 
tivated goodness of fit tests in the spike train analysis field had to wait until the 
beginning of this century with a paper by [3, Brown et al, 2002]. This latter work 
introduced in neuroscience one of the tests proposed by [ , Ogata, 1988]. All these 
tests are based on the fundamental "time transformation" [16] / "time rescaling" [3] 
result stating that is the conditional intensity model is correct, then upon mapping 
the spike times onto their integrated conditional intensity a homogeneous Poisson 
process with rate 1 is obtained. The tests are then all looking for deviations from 
homogeneous Poisson process properties. 

We will not address the estimation of the conditional intensity in this paper in or- 
der to keep its length within bounds. Our own take at this problem will be presented 
elsewhere (Pouzat, Chaffiol and Gu, in preparation). Simple examples of conditional 
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intensities will be used for simulations and applied to actual data. We will focus in- 
stead on a new goodness of fit test based on a straightforward application of Donsker's 
Theorem [ ]. We will show that the test has good finite sample properties and sug- 
gest that it is a useful complement to Ogata's tests. We will moreover argue based on 
real data examples in favor a systematic use of several tests. 

The paper is divided as follows. In Sec. 2, an elementary justification of the time 
transformation is presented. This proof is obtained as a limit of a discrete problem 
arising from a uniform binning of the time axis. This discretized problem is more- 
over the canonical setting appearing in (nearly) every presently proposed conditional 
intensity estimation method. A simple illustration of the time transformation is pre- 
sented using simulated data. In Sec. 3, the "Ogata's tests battery" is presented, dis- 
cussed and illustrated with the simulated data set of the previous section. In Sec. 4, 
a new test, the "Wiener process test", is proposed based on Donsker's Theorem. The 
Theorem is first stated and a description of how to use it on time transformed data is 
given. A "geometrical" interpretation is sketched using the simulated data set. The 
finite sample properties of the test are studied next with Monte Carlo simulations. 
In Sec. 5, real data are used to illustrate the performances of the Ogata's tests and of 
the Wiener process test. Sec. 6 summarizes our findings and discusses them. 

2 An heuristic approach to the time transformation 

A simple and elegant proof of the time transformation / time rescaling theorem ap- 
pears in [3]. We are going to present here an heuristic approach to this key result 
introducing thereby our notations. Following [2] we define three quantities: 

• Counting Process: For points {tj} randomly scattered along a line, the count- 
ing process N(t) gives the number of points observed in the interval (0, t]: 

N(t) = D-fo with < tj < t} (1) 

where fj stands for the number of elements of a set. 

• History: The history, Tit, consists of the variates determined up to and includ- 
ing time t that are necessary to describe the evolution of the counting process. 
Tit can include all or part of the neuron's discharge up to t but also the discharge 
sequences of other neurons recorded simultaneously, the elapsed time since the 
onset of a stimulus, the nature of the stimulus, etc. One of the major problems 
facing the neuroscientist analysing spike trains is the determination of what 
constitutes Tit for the data at hand. A pre-requisite for practical applications of 
the approach described in this paper is that Tit involves only a finite (but possi- 
bly random) time period prior to t. In the sequel we will by convention set the 
origin of time at the actual experimental time at which Tit becomes defined. For 
instance if we work with a model requiring the times of the last two previous 
spikes, then the t = of our counting process definition above is the time of the 
second event / spike. If we consider different models with different Tit, the time 
origin is the largest of the "real times" giving the origin of each model. 

• Conditional Intensity: For the process N and history Tit, the conditional in- 
tensity at time t is defined by: 

»(*!«,) = lim P"*WM + q-Ar W = i|K,} 

510 5 
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As far as we know, most of the presently published attempts to estimate the con- 
ditional intensity function involve a discretization (binning) of the time axis [10,2,5, 
11,23, 17,22]. We will follow this path here and find the probability density of the 
interval between two successive events, Ij = tj+i — tj. Defining 5 = ti+1 ~ tj , where 
K e W. We first write the probability of the interval as the following product: 

Px{I j =t j+x -t j ) = Pr{N(t j + S)-N(t j ) = 0\H tj }- 

• Pr{N(tj + 25)- N(tj + S)=Q\ H tj+S } ■ ■ ■ 

■ ■ • Pr{N(t, + K 5) - N(tj + (K - 1) 5) = | H t ^ K _ X) s } ■ 

■ PT{N(tj + (K+ 1)5) - N(tj + K5) = 1 | H tj+KS } (3) 

Then using Eq. 2 we can, for K large enough, consider only two possible outcomes 
per bin, there is either no event or one event. In other words we interpret Eq. 2 as 
meaning: 

Prob{iV (t, t + 5) - N(t) = | H t } = 1 - A(i | H t ) 6 + o(5) 
Proh{N{t,t + 5) - N(t) = 1 | H t } = X(t | H t ) 5 + o(5) 
Pvoh{N (t, t + S) - N(t) >l\H t } = o(S) 

where o(<5) is such that lim^^o ^ = 0. The interval's probability (Eq. 3) becomes 
therefore the outcome of a sequence of Bernoulli trials, each with an inhomogeneous 
success probability given by Aj 5 + o(5), where, \ = X(tj + i 5 \ Htj+is) and we get: 

K 

Pv{Ij = t j+1 - tj} = ( J](l - \ k 5 + o(S))) {X K +i S + o(5)) (4) 

k=l 

We can rewrite the first term on the right hand side as: 

K K 



Y[(l-\ k 5 + o(5)) = explogJ](l-A fc «5 + o(5)) 

k=l 

K 

exp ^2 log(l - X k 5 + o(6)) 

k=l 
K 

exp^(-A fc 5 + o(S)) 

k=l 

K 

exp(- ^XkS) ■ exp (K o(5)) 



k=l k=l 

K 



k=l 



Using the continuity of the exponential function, the definition of the Riemann's in- 
tegral, the definition of 5 and the property of the o() function we can take the limit 
when K goes to oo on both sides of Eq. 5 to get: 

lim TT (1 - A fc 5 + o(<5)) = exp - / ' X(t \ H t ) dt (5) 

k=l Jt o 

And the probability density of the interval becomes: 



Pv{Ij=t J + 1 -tj} 



lim J t JX = A (^+i I Ht j+1 ) exp - / X(t\ H t ) dt (6) 



K 



■j+i 
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If we now define the integrated conditional intensity function by: 



A(t) = f A( 

Ju=0 



u I H u ) du (7) 

We see that A is increasing since by definition (Eq. 2) A > 1 . We see then that the 
mapping: 

t g R + * -> A(t) G M+* (8) 

is one to one and we can transform our {t±,..., t n } into {Ai = A(ti), . . . , A n = A(t n )}. 
If we now consider the probability density of the intervals tj+i — tj and Aj + i — Aj we 
get: 

rtj+i 

p(t j+ i - tj) dt j+1 = \(tj+i | H tj+1 ) exp ( - / \(t | H t ) dt) dt j+1 

= ^l^^+iexp-lA^+O-A^)) 

= dA j+1 exp - (A j+ i - Aj) (9) 

That is, the mapped intervals, Aj + \ — Aj follow an exponential distribution with rate 
1. This is the substance of the time transformation of [16] and of the time rescaling 
theorem of [ ]. This is also the justification of the point process simulation method 
of [9, Sec. 2.3, pp 281-282]. 

An illustration 

We illustrate here the time transformation with simulated data. This is also the occa- 
sion to give an example of what an abstract notation like, \(tj \ Ht d ), translates into 
in an actual setting. We will therefore consider a counting process whose conditional 
intensity function is the product of an inverse-Gaussian hazard function and of the 
exponential of a scaled x 2 density. That would correspond to a neuron whose spon- 
taneous discharge would be a renewal process with a inverse-Gaussian inter spike 
interval (isi) distribution and that would be excited by "multiplicative" stimulus (this 
type of Cox-like model was to our knowledge first considered in neuroscience by [10]). 
More explicitly our conditional intensity function is defined by: 

h, a(x) = -lip- (12) 
1 - F IG (x) 

s(t) = pf xl (m(t-t )) (13) 

A(t | H t ) = h IG (t - t t ) exp (s{t)) (14) 

where, is the cumulative distribution function of a standard normal random vari- 
able and / 2 stands for the probability density function of a x 2 random variable with 5 
degrees of freedom and t\ is the occurrence time of the preceding spike. In this rather 



1 Actual neurons exhibit a "refractory period", that is, a minimal duration between two successive 
spikes. One could therefore be tempted to allow A = on a small time interval following a spike. We 
can cope with this potential problem by making A very small but non null leading to models which 
would be indistinguishable from models with an absolute refractory period when applied to real world 
(finite) data. 
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simple case the history is limited to: Tit = max{t, : tj < t}. The graph of X(t \ Tit) 
corresponding to 10 s of simulated data is shown on Fig. 1 A, together with the data. 
The discontinuous nature of the conditional intensity function appears clearly and 
is a general feature of non-Poisson counting processes. Fig. 1 B, shows the graph of 
the counting process together with the integrated conditional intensity function, A(t). 
Fig. 1 C, illustrates the time transformation per se, the counting process is now de- 
fined on the "A scale" and the integrated conditional intensity is now a straight line 
with slope one on this graph. 

3 The Ogata's tests 

As soon as we have a model of X(t \ Tit) that we think (or hope) could be correct 
we can, following Ogata [16], exploit the time transformation of the previous sec- 
tion to generate goodness of fit tests. We start by mapping {h, . . . ,t n } onto {Ai = 
A(ti), . . . , A n = A(t n )}. Most of the time this mapping requires a numerical integra- 
tion of X(t | Tit). Then if our model is correct, the mapped intervals, A J+ i - Aj, should 
be iid from an exponential distribution with rate 1 and (equivalently) the counting 
process N(A) should be the realization of a homogeneous Poisson process with rate 1. 
In [16], Ogata introduced the five following tests: 

1. If a homogeneous Poisson process with rate 1 is observed until its nth event, 
then the event times, {Aj}™^ 1 , have a uniform distribution on (0, A n ) [6, chap. 
2]. This uniformity can be tested with a Kolmogorov test. This is the first Ogata 
test [ , Fig. 9, p 19]. It's application to our simulated data is shown on Fig. 2 
A. In the sequel we will refer to this test as Ogata's first test or as the "uniform 
test". 

2. The u k defined, for k > 1, by: 

u k = 1 - exp ( - (A fc - A fe _i)) 

should be iid with a uniform distribution on (0, 1). The empirical cumulative 
distribution function of the sorted {u/J can be compared to the cumulative dis- 
tribution function of the null hypothesis with a Kolmogorov test. This test is 
attributed to Berman in [ ] and is one of the tests proposed and used by [3]. 
It's application to our simulated data is shown on Fig. 2 B. This tests implies 
a sorting of the transformed data which would destroy any trace of serial cor- 
relation in the u k s if one was there. In other words this tests that the -u^s are 
identically distributed not that they are independent. This latter hypothesis has 
to be checked separately which leads us to the third Ogata test. 

3. A plot of u k+1 vs Uk exhibiting a pattern would be inconsistent with the homo- 
geneous Poisson process hypothesis. This graphical test is illustrated in Fig. 2 
C for our simulated data. A shortcoming of this test is that it is only graph- 
ical and that it requires a fair number of events to be meaningful. In [20] a 
nonparametric and more quantitative version of this test was proposed. 

4. Ogata's fourth test is based on the empirical survivor function obtained from the 
tifcS. It's main difference with the second test is that only pointwise confidence 
intervals can be obtained. That apart it tests the same thing, namely that the 
marginal distribution of the u k s is exponential with rate 1. Like the second 
test it requires sorting and will fail to show deviations from the independence 
hypothesis. Since this test being redundant with test 2 it is not shown on Fig. 2. 



5 



Conditional intensity and events sequence 




2 4 6 8 10 50 100 150 

Time (s) A 

Figure 1: Time transformation illustration. From simulated data with the thinning 
method [15] according to the conditional intensity function defined by Eq. 14. Pa- 
rameters used n = 0.075 and a 2 = 3 in Eq. 10 and 11; to = 4, m = 5 and p = 20 in 
Eq. 13. A, Graph of the conditional intensity function together with the spike train 
represented by a sequence of small vertical bars (one for each spike at each spike 
time) at the bottom of the graph. The conditional intensity function is discontinuous, 
exhibiting a jump at each spike occurrence. It is left continuous and converges to 
on the right-hand side. B, the integrated conditional intensity function, A (continuous 
curve) and the counting process, N (step function) as a function of time. C, Same as 
B after time transformation. 
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Uniform on A Test Berman's Test 





Figure 2: Four of the five Ogata's tests applied to the simulated data of Fig. 1. 



5. The fifth test is obtained by splitting the transformed time axis into K w non- 
overlapping windows of the same size w, counting the number of events in each 
window and getting a mean count N w and a variance V w computed over the K w 
windows. Using a set of increasing window sizes: {w%, . . . ,wl} a graph of V w 
as a function of N w is build. If the Poisson process with rate 1 hypothesis is 
correct the result should fall on a straight line going through the origin with 
a unit slope. Pointwise confidence intervals can be obtained using the normal 
approximation of a Poisson distribution as shown on Fig. 2 C. 



4 A new test based on Donsker's Theorem: The Wiener 
Process Test 

As we have seen in the previous section each Ogata's test considered separately has 
some drawback. This is why Ogota used all of them in his paper [16]. It seems 
therefore reasonable to try to develop new tests which could replace part of Ogata's 
tests and / or behave better than them in some conditions, like when the sample size 
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is small. We are going to apply directly Donsker's Theorem to our time transformed 
spike trains as an attempt to fulfill this goal. 

Following Billingsley [1, p 121], we start by defining a functional space: 

Definition of D Let D = D[0, 1] be the space of real functions x on [0, 1] that are 
right-continuous and have left-hand limits: 

(i) For < t < 1, x(t+) = linisji x(s) exists and x(t+) = x(t). 

(ii) For < t < 1, x(t-) = lim^ x(s) exists. 

Donsker's Theorem Given iid random variables £i, £2, • • • , £n, ■ • • with mean and 
variance a 2 defined on a probability space (fi, F, P). We can associate with these 

random variables the partial sums S n = £1 H 1- £ n - Let X n {uj) be the function in D 

with value: 

X t » = ^ 7 =5 LtrfJ H (15) 

at t, where ["^J is the largest integer < nt. Donsker's Theorem states that: The 
random functions X n converge weakly to W, the Wiener process on [0, 1]. 
See [ , pp 146-147] for a proof of this theorem. 

Using Donsker's Theorem If our model is correct then the A J+ i — Aj are iid ran- 
dom variables from an exponential distribution with rate 1. Since the mean and 
variance of this distribution are both equal to 1, if we define: 

& = A j+1 -Aj-1 (16) 

the £jS should be iid with mean and variance 1. Once we have observed the AjS and 
therefore the £jS we can built the realization of Eq. 15. The latter should "look like" 
a Wiener process on [0, 1]. Notice that the realization of Eq. 15 is nearly the same 
as the function obtained by subtracting N(A) from A on Fig. 1 C, before dividing the 
abscissa by n and the ordinate by y/n as illustrated on Fig. 3. This is also seen by 
looking at the expression of the partial sums in this particular case: 

S n = A n+1 - Ai - n (17) 

The next element we need in order to get a proper test is a way to define a tight 
region of [0, l]xR where, say, 95% of the realizations of a Wiener process should be. 
This turns out to be a non trivial but luckily solved problem. Kendall et al [12] have 
indeed shown that the boundaries of the tightest region containing a given fraction 
of the sample paths of a Wiener process is are given by a so called Lambert W func- 
tion. They have moreover shown that this function can be very well approximated 
by a square root function plus an offset: a + b \/t. We next need to find the a a and 
b a such that the realizations of a Wiener process are entirely within the boundaries 
with a probability 1 — a. There is no analytical solution to this problem but efficient 
numerical solutions are available. Loader and Deely [13] have shown that the first 
passage time of a Wiener process through an arbitrary boundary c(t) is the solution of 
a Volterra integral equation of the first kind. They have proposed a "mid-point" algo- 
rithm to find the first passage time distribution. Their method does moreover provide 
error bounds when c(t) = a + b \ft. Then using the "symmetry" of the Wiener process 
with respect to the time axis we have to find a a and b a such that the probability to 
have a first passage < 1 is 
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Figure 3: Illustration of the "mapping" to a Wiener process. A, The counting process 
has been subtracted from the integrated intensity of Fig. 1 C, on the transformed 
time scale. B, The corresponding realization of the X t n function of Eq. 15. 



Wiener Process Test 



CM 
I 



CO 
I 




0.2 



0.4 



0.6 



0.8 



1.0 



Figure 4: Illustration of the "Wiener process test". The same graph as Fig. 3 B is 
shown with 95% and 99% boundaries appearing as dotted curves. The boundaries 
have the form: a + b \ft, the coefficients are given in the text. 
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Using an integration step of 0.001 and, a . 05 = 0.299944595870772, 6 .05 = 2.34797018726827, 
we get: 0.9499 < P .95 < 0.9501. Using o .oi = 0.313071417065285, 6 ,oi = 2.88963206734397, 
we get: 0.98998 < Po.99 < 0.99002. In the sequel we will refer to the comparison of the 
sample path of X™ in Eq. 15 with boundaries whose general form is: a a + b a y/t as a 
Wiener process test. 

Finite sample properties We now estimate the empirical coverage probability of 
our two regions for a range of "realistic" sample sizes, from 10 to 900, using Monte 
Carlo simulations. The results from 10000 simulated "experiments" are shown on 
Fig. 5 A. The empirical coverage probability of the region with a nominal 95% confi- 
dence is at the nominal level across the whole range studied. The empirical coverage 
probability of the region with a nominal 99% confidence does slightly worse. The em- 
pirical coverage is roughly 98% for a sample size smaller than 100, 98.5% between 
100 and 300 and reaches the nominal value for larger sample sizes. The overall per- 
formances of the test are surprisingly good given the simplicity of its implementation. 

Since we are also interested in carrying out multiple testing using for instance 
Ogata's first test, Berman's test and our "Wiener process test", we studied the depen- 
dence of false negatives under the three possible pairs formed from these three tests. 
The results are shown on Fig. 5 B. Here we see that Berman's and Ogata's first tests 
are independent {i.e., the empirical results are compatible with the hypothesis of in- 
dependence). But neither of these tests is independent of the "Wiener process test" 
when a 95% confidence level (on each individual test) is used. In both cases (Berman 
and Ogata's first test) when a trial is rejected at the 95% level it is less likely to be 
also rejected by the "Wiener process test". If we combine the three tests at the 99% 
level and decide to reject the null hypothesis if any one of them is rejected, we obtain 
an empirical coverage probability of ~ 0.96 for a sample size between 10 and 100. 
Above this size the probability becomes ~ 0.97. 

A remark It is clear from the hypothesis of Donsker's theorem that the "Wiener 
process test" cannot replace Ogata's second test (Berman's test). As long as the dis- 
tribution of the mapped intervals has the same first and second moments than an 
exponential distribution with rate 1, the data will pass the test. 

5 Real data examples 

We illustrate in this section the performances of the tests on real data. The data, 
like the software implementing the test and the whole analysis presented in this 
article are part of our STAR (Spike Train Analysis with R) package for the R 2 software. 
The data were recorded extracellularly in vivo from the first olfactory relay of an 
insect, the cockroach, Periplaneta americana. Details about these recordings and 
associated analysis methods can be found in [ ,20]. The names of the data sets used 
in this section are the names used for the corresponding data in STAR. The reader 
can therefore easily reproduce the analysis presented here. We do moreover provide 
with STAR a meta-file [21,7], called a vignette in R terminology, allowing the reader to 
re-run exactly the simulation / data analysis of this article. 

5.1 Spontaneous activity data 

Two examples of spike trains recorded during a period of 60 s of spontaneous activity 
are shown on Fig. 6. The spike trains appear as both a counting process sample path 

2 http : / / www .r-project.org 
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Figure 5: A, Empirical coverage probabilities of the 95 and 99% confidence regions 
for the Wiener test. Each dot on the graph was obtained from 10000 simulated "ex- 
periments". Each experiment was made of a number, n, of draws given on the ab- 
scissa. Each draw was the realization of an exponential distribution with rate 1. The 
procedure described by applying Eq. 16 before Eq. 15 was then applied. The result- 
ing step function was compared, at the steps locations with the boundaries values. 
The grey band correspond to 95% confidence bands for a binomial distribution with 
10000 draws and a success probability of 0.95 (bottom) and 0.99 (top). B, Number 
of joint rejections of paired tests. Using the same simulations as in A, the number 
of joint rejection of the "Wiener process test" and of the "Berman's test" (circles), of 
the "Wiener process test" and of the "Uniform test" (i.e., Ogata's first test, triangles), 
of the "Berman's test" and of the "Uniform test" (crosses) at the 95% level (dotted 
lines) and at the 99% level (continuous line). Gray areas, 95% confidence region for a 
binomial distribution with 10000 trials and a rejection probability of 0.05 2 (top) and 
0.01 2 (bottom). The p values of the Kolmogorov statistic of the "Berman" and "Uni- 
form" tests were obtained with the exact method of [14] and implemented in function 
ks . test of R. 
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Neuron 3, data set: e060517spont 



Neuron 1, data set: e060824spont 



o 




10 20 30 40 50 60 10 20 30 40 50 



Time (s) Time (s) 

Figure 6: Two actual spike trains. Each plot shows the counting process sample path 
(N(t)) associated with a train. The spikes density and their number are too high for 
the discontinuities of N(t) to systematically appear. A "raster plot" is drawn at the 
bottom of each plot. A raster plot is a spike train representation where the occurrence 
time of each spike is marked by a tick (see [ ] for details). 

and as a "raster plot" on this figure. A renewal model was fitted with the maximum 
likelihood method to each train. The data from neuron 3 of data set e060517spont 
were fitted with an inverse-Gaussian model, while the ones of neuron 1 of data set 
e060824spont were fitted with a log logistic model. The spike times where trans- 
formed (Eq. 7) before applying Ogata's test 1, 2, 3, 5 and the new Wiener test. The 
results are shown on Fig. 7. 

The time transformed spike train of neuron 3 of data set e060517spont passes 
Berman's test (Fig. 7 A2) as well as Ogata's third test (Fig. 7 A3) but fails to pass 
Ogata's tests 1 and 5 (Fig. 7 Al and A4) as well as the new Wiener process test (Fig. 7 
A5). The time transformed spike train of neuron 1 of data set eO 60824 spont passes 
Ogata's first test (Fig. 7 Bl) as well as Berman's test (Fig. 7 B2), but fails Ogata's 
tests 3 (Fig. 7 B3) and 5 (Fig. 7 B4) as well as the new Wiener process test (Fig. 7 B5). 



5.2 Stimulus evoked data 

The odor evoked responses of neuron 1 of data set e070528citronellal are illus- 
trated here. This is the occasion to use a more sophisticated intensity model, albeit 
still a wrong one (see bellow). The "raw data" are shown on Fig. 8. An inhomogeneous 
Poisson model was fitted to trial 2 to 15 using a smoothing spline approach [8, 20]. 
This estimated time dependent intensity was then used to transform the spike times 
of the first trial. The tests are shown on Fig. 9. In this case the train contains 98 
spikes making Ogata's test 5 barely applicable (not shown on figure). Ogata's test 3 
(Fig. 9 C) is hard to interpret since with few spikes and a sparse filling of the graph, 
the presence or absence of a pattern is not obvious. Ogata's test 1 (Fig. 9 A) and 
Berman's test (Fig. 9 B) are passed but the Wiener process test (Fig. 9 D) is not. 
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Uniform on A Test 
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Berman's Test 





B2 



Berman's Test 




A5 



Window Length 
Wiener Process Test 



B5 



Window Length 
Wiener Process Test 




Figure 7: The tests applied after time transformation. A, neuron 3 from data set 
eO 60517 spont fitted with an inverse-Gaussian renewal model. B, neuron 1 from 
data set e060824spont fitted with a log logistic renewal model. Notice that the 
tests are here applied directly to the data used to fit the model. The domains defined 
by the boundaries (dotted lines) have therefore not exactly their nominal coverage 
value. 
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Figure 8: An example of odor evoked responses. Neuron 1 of data set 
e070528citronellal is used. 15 citronellal odor puffs were applied successively. 
Each puff was 0.5 long [4, 20] (the gray area corresponds to the opening time of the 
valve delivering the odor). Odor presentations were performed 1 mn apart. The spike 
trains of each trial are represented as raster plots. The first trial is at the bottom of 
the graph, the 15th on top. 
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6 Conclusions 



Ogata [ ] introduced a tests battery for counting process models. He made moreover 
rather clear through his use of the tests that more than a single one of them should be 
used in practice. Brown et al [3] introduced in the neuroscience literature Berman's 
test (Ogata's second test) which has become the only test used (when any is used) 
in this field. We could also remark that only one case of a model fitted to real data 
actually passing this test has ever been published [3, Fig. 1A]. Clearly our examples 
of Fig. 7 A2 and B2 and of Fig. 9 B, would multiply this number by four. But we 
have to stick to Ogata's implicit message. Our choice of examples reiterates this 
point showing cases where Berman's test is passed. We do not want to imply that 
Berman's test is a "bad" test but to warn the neuroscience community that misleading 
conclusions can be obtained when only this test is used. 

We have also introduced the "Wiener process test" a direct application of Donsker's 
Theorem to the time transformed spike trains. The test is simple to implement 
(Eq. 16 and 15) and exhibits good finite sample properties (Fig. 5). It is clearly not 
a replacement for the Ogata's tests battery but a complement (Fig. 7 and 9). It's fi- 
nite sample properties make it particularly useful in situations were a small number 
of spikes makes the interpretation of Ogata's tests 3 and 5 potentially ambiguous 
(Fig. 9). 
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